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We study the Casimir interaction between a metallic cylindrical wire and a metallic spherical 
particle by employing the scattering formalism. At large separations, we derive the asymptotic form 
of the interaction. In addition, we find the interaction between a metallic wire and an isotropic 
atom, both in the non-retarded and retarded limits. We identify the conditions under which the 
asymptotic Casimir interaction does not depend on the material properties of the metallic wire and 
the particle. Moreover, we compute the exact Casimir interaction between the particle and the wire 
numerically. We show that there is a complete agreement between the numerics and the asymptotic 
energies at large separations. For short separations, our numerical results show good agreement 
with the proximity force approximation. 



I. INTRODUCTION 

Casimir forces contribute significantly to the effective 
interaction of micro- and nanometer sized structures [1] . 
For identical objects or mirror symmetric configurations, 
this type of interaction is attractive [2] and can cause 
stiction in micro-motors and other similar structures [3]. 
More generally, if the permittivities of the objects are 
higher or lower than those of the surrounding medium, 
any equilibrium position of the objects is unstable due 
to the Casimir interactions [4]. Therefore, a good quan- 
titative understanding of such forces is a key parameter 
in the design and manufacturing of micro-mechanical de- 
vices. 

It is important to study the Casimir forces for different 
shapes as they strongly depend on the geometry and ma- 
terial properties [TJ O H] . Technically, investigating the 
interplay between the shape and material effects is quite 
involved. The scattering formalism provides a powerful 
tool to calculate the Casimir interaction between objects 
of general shape and material properties [7 . There is 
much recent research activity based on the scattering for- 
malism, e.g., for edges and tips [HUH]) anisotropic parti- 
cles [TU] , wires and plates [PflfB] , and spheres and plates 

An important geometry which has not yet been in- 
vestigated in detail, consists of a wire and a particle 
(atomic or macroscopic). In the plane-particle geome- 
try this force is known as Casimir-Polder (CP) interac- 
tion [IB]. Our study of the wire-particle case is moti- 
vated by theories [IT] and experiments jTSJ IH] on the 
two-dimensional quantum scattering of neutral atoms or 
molecules at wires or nanotubes. In an early work, the 
interaction between a filament and an isotropic atom has 
been studied for perfectly and non-perfectly conducting 
metals [20]. Later, Eberlein and Zietal studied the inter- 
action between a neutral atom and a perfect metal cylin- 
der, using perturbation theory [3TJ [22]. Recently, the 



Casimir energy for a polarizable micro-particle and an 
ideal metal cylindrical shell has been computed using the 
Green's function technique [23j . The focus of previous 
studies were mainly on the interaction between a metal 
wire and a perfect metal particle or an atom. There- 
fore, the influence of material properties of the spherical 
particles on the energy remains to be studied in detail. 

In this work, we study the Casimir interaction between 
a metallic spherical particle and cylindrical metallic wire 
where the latter is described by the Drude, plasma or 
perfect metal model. Using the scattering formalism, 
we derive a general expression for the Casimir interac- 
tion between the particle and cylinder. From this gen- 
eral expression we determine the behavior of the Casimir 
energy in various limiting cases (separation regimes) an- 
alytically, and numerically over a wide range of separa- 
tions. Interestingly, we find ranges of distances in which 
the Casimir interaction does not depend on the mate- 
rial properties of the metallic wire. In contrast, we find 
that the interaction depends in general on the mate- 
rial properties of the metallic particle at all separations. 
An exception is the plasma sphere with a plasma wave- 
length that is much smaller than the size of the sphere for 
which the Casimir interaction is universal at asymptoti- 
cally large distances. At short separations, we compute 
the exact Casimir interaction numerically and compare 
it both with the asymptotic results and the prediction of 
proximity force approximation (PFA). In both limits we 
obtain good agreement. 

The structure of the rest of the paper is as follows: In 
Sec. [IT] we review the scattering approach and derive the 
elements that are needed for computing the interaction 
between a wire and a particle. In Sec. |III| the large- 
distance asymptotic interaction between a metallic wire 
and a particle (metal sphere and isotropic atom) is de- 
rived for the perfect metal, plasma and Drude models. 
In Sec. |IV[ the exact Casimir interaction is computed 
numerically and compared with the asymptotics expan- 



2 



sions. Section [V] is dedicated to the interaction at short 
separations where the PFA is expected to become reli- 
able. 



II. METHOD 

We consider a cylindrical wire separated by a distance 
d from a spherical particle. We use the scattering formal- 
ism to calculate the Casimir energy between the cylinder 
and sphere [7j. In general, the Casimir free energy be- 
tween two objects at the temperature T is given by 



In det [1 



(1) 



n=0 



where 1 is the identity matrix, K n is the Matsubara wave 
number n n = 2nnkBT/hc and the matrix N factorizes 
into the scattering amplitudes (T-matrices) as well as 
translation matrices which describe the coupling between 
the multipoles on distinct objects. The primed sum de- 
notes that the contribution of n = has to be weighted 
by a factor of 1/2. 

At zero temperature, the primed sum in Eq. fl) is re- 
placed by an integral along the imaginary frequency axis, 



he f°° 
£ = — / d/clndet(l 



N) 



(2) 



with k the Wick-rotated frequency. The elements of the 
matrix N for electric (E) and magnetic (M) polarizations 
(a, /3, 7 = E, M) and cylindrical wave functions m and 



m are 



N afj \ " \ "* T a i 

k z m,k'^m' / j / t s,k z m,k' z m" 

*y=E,M m" — — oo 



where k z is the wave number along the z-axis, T c and T s 
are the T-matrices of the cylinder and sphere in cylindri- 
cal basis, respectively. The translation matrix U sc relates 
regular cylindrical vector waves to outgoing ones. 

The translation matrices do not couple different po- 
larizations and for both E and M polarizations, their 
matrix elements are given by 



ZC, = (-1)" K n _ n ,(pd), 

7 yes / -1 \n— n 7 jsc 

U nn> — \ l ) U nn> j 



(4) 



where p = \f k 2 + k% and K n (x) is the modified Bessel 
function of the second kind. Note that the N-matrix ele- 
ments in Eq. ([3 ) are written in cylindrical basis to avoid 
the complicated form of the translation matrices in the 
spherical basis [7J. 



The T-matrix of the sphere in cylindrical basis is de- 
rived in Appendix [B] and is given by 



s,k z m,k'm 



£=max(l,|m|) f) 



xD 
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(5) 



where L is the length of the cylinder, I is the quan- 
tum number of the spherical electromagnetic waves, /3 = 
E, M is the electromagnetic polarization and D is the 
conversion matrix from the cylindrical to spherical basis. 
The elements of the conversion matrix are given in Ap- 
pendix [cj. Note that in Eq. (j5j, T^ tra is the T-matrix of 
the sphere in the spherical basis. 

To obtain the Casimir energy from Eq. we plug 
Eqs. Q and Q into Q and use the identity det(l - 
AB) = det(l — BA). The N-matrix for the energy be- 
tween the sphere and the cylinder is rewritten as 



1 r c 
_ 1 V T a \ 

•"»•''»'' 47T 2 K s ' tm J 



xTl^U^a-^y). ( 6 ) 



with 



~ 77 ' _ X ■> j. cs rr^ff' j JSC /y\ 
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In this work, to study the impact of the material prop- 
erties of the metallic objects on the Casimir interaction, 
we employ the plasma, Drude and perfect metal dielec- 
tric properties with the constant magnetic permeability 
/i = 1. The Drude model dielectric response is given by 



e{icK, A p , A CT ) = 1 + 



(2tt) s 



(A pK ) 2 + X aK /2 



(8) 



where A p is the plasma wave length and A CT = 2irc/a is the 
length scale associated with the conductivity a. Equa- 
tion Q , reproduces the plasma model dielectric function 
for \fj — 0. Note that the material properties of the 
sphere and the cylinder enter to the calculations through 
the T-matrices, see Appendix [A] 



III. LARGE-SEPARATION REGIME: 
ASYMPTOTIC CASIMIR ENERGY 

In this section, we study the large separation asymp- 
totic behavior of the Casimir interaction between a par- 
ticle and a wire. We consider a spherical particle with 
radius R s and a cylindrical wire with radius R c . In order 
to find the large-distance (d ^> R C ,R S ) asymptotic form 
of the Casimir interaction, one has to find the behavior 
of the T-matrices in the low- frequency limit. 
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A. Asymptotic Behavior of T-matrices 

1. T -matrix of a wire 

In this part, we obtain the asymptotic form of the T- 
matrix elements of a wire at large separations. Using the 
dielectric function given in Eq. (Jsj), we find the T-matrix 
element of the wire for E polarization and n = at small 
frequencies (k <C 1, k z / k fixed) is 



P 



C(k)- P 2 \n{pR c 



(9) 



where p = K 2 + k 2 . The parameter C(k) depends on 
the dielectric properties of the wire. For a perfect metal 
wire C(k) = and for a plasma wire with the plasma 
wavelength A p , C{n) w A 2 k 2 / '(2ir 2 R 2 ) if the plasmon os- 
cillations cannot build up transverse to the wire axis as 
the diameter is too small, i.e., R c <C A p . In the op- 
posite limit we approximately reproduce the T-matrix 
of a perfect metal wire, i.e. C(k) 0. For a Drude 
wire with the conductivity a and the characteristic length 
A CT , C(k) = A ff K/(47r 2 i?2) if k < A^/A 2 , 1/A CT . The first 
condition (n <C A CT /A 2 ) means that the Drude behavior 
dominates over the plasma behavior, equivalent to the 
fact that in the dominator of Eq. pL the first term is 
much smaller than the second one. The second condition 
(k <C 1/A ct ) ensures that the Drude dielectric function is 
much larger than one, i.e., the metallic behavior is dom- 
inant PUDl]. 

At large separations, T E E elements dominate over the 



other T-matrix elements since T, 



EM 



rpM M 
1 



rpME 

2 



and for n ^ partial waves T n 



= 0, 

K 2\n\ 



Note that for Drude cylinders T^ E 
and perfect metal cylinders T EE ~ 



k while for plasma 



1. 



2. T-matrix of a particle 

The T-matrix elements of a spherical particle have 
a different scaling compared to the ones for the cylin- 
drical wire. For the plasma and perfect metal spheres 
T Sl i m ~ K 2i+1 . For Drude spheres the T-matrix elements 
for the electric E and magnetic M polarizations scale 
differently: T s A j m - k 2 ( £+1 ) and T s E im - n 2l+1 . There- 
fore, the asymptotic behavior of the T-matrix at large 
separations is dominated by the I = 1 elements. The 
asymptotic form of the T-matrix elements for the mag- 
netic polarization M and £ = 1 depends on the mate- 
rial properties of the sphere. While for a perfect metal 
sphere, we have 



s,lm 



k 3 R 3 J3 , 



(10) 



for a plasma sphere with the plasma wavelength Ap, 
obtain 



we 



1 I i (2nR s /X' p ) 
'3/i(27riyA' p ) 



K 3 Rt 



(11) 



Note that in the limit of perfect conductivity A p — > 0, we 
reproduce Eq. (10). 



For a Drude sphere with the conductivity a' and the 
characteristic length A^, 



rpM ^ _ & 71 " 2 4 R 4 

s ' lm ~ 45 A' 



(12) 



However, the asymptotic form of the T-matrix elements 
for E polarization and I — 1, up to the leading order, 
does not depend on the material properties, and is the 
same for the perfect metal, plasma and Drude models, 



T E at -k 3 P 3 
J s,lm ~ g ll s 



(13) 



Note that for the Drude sphere T s A { m ~ k 4 , and thus we 
add the sub-leading term to the expansion of T E lm in 
Eq. ( 13 ) and obtain 



rpE 

J s,lr, 



3 D 3 



-K 6 R 



4tt 2 JL K Hs ' 
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where the sub-leading term contains the material prop- 
erties of the Drude wire. 



3. T-matrix of an Atom 

The above approach can also be used to calculate the 
Casimir energy between an atom and a wire. To this 
end, we consider a neutral two-level atom in the ground 
state, with the transition frequency wio [16] . We assume 
that the distance from the atom to the wire d is much 
larger than the radius of the wire R c , i.e. d/R c 1. 
Moreover, we assume that the atom is isotropic and does 
not have magnetic polarizability. In the isotropic-dipole 
approximation, the only nonzero element of the T-matrix 
reads 



Ti 



r 



(k) k 3 



where a E , the electric polarizability, is given by 

on 



<x e {k) 



1 + K 2 d 2 ' 



(15) 



(16) 



with dw = c/wio, «o = fw^ /(mu}f ), e the electron 
charge, m the mass and /io the oscillator strength of the 
— >• 1 transition. 



B. Asymptotic Energy Expression 

In this subsection, using the asymptotic T-matrix ex- 
pressions and Eq. ([6]), we derive the Casimir energy at 
large separations. Considering lndet = Trln, we expand 
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the integrand in Eq. ^ in powers of N for re <C 1 and 
k z /n fixed and find 



£ 



he 
"2^ 



dn TrN . 



As discussed above, in the limit d ^> R S ,R C only the 
I = 1 terms contributes to the sphere T-matrix, T s ^ m . 
Therefore, only the partial wave numbers m — 0, ±1 have 
to be taken into account to obtain the matrix N. Using 
Eq. (|6|, we find 



TrN 



7 m— — 1 



r 2 M . 



47T 2 K 
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V T 7 

/ j s,lm 



7 m—— 1 



w/v^ i-^im^ ,k z mE mm 



xDlmE,lmj^-2Sj, E ), (18) 



with D the modified conversion matrix with real ele- 
ments. The modified conversion matrix is related to the 
original one by B em . k;cm = (-1 ) t ~ m {-iY +m ~ 1 'Di ni> k % m 
, see Appendix [c] Furthermore, the matrix T^m m 
Eq. ( 18 ) is equal to 



-\EE 



c.fcjO - 



(19) 



Note that in Eq. (19 1, T^f ±1 are neglected as they scale 
with higher powers of re. 



Inserting Eq.(19) into Eq. (18) and performing the 
sums, we find 
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The modified conversion matrix elements in Eq. ( 20 1 are 
(see Appendix [C] for all details ), 

D 10E ,k a 0M = 0, D WM .k,0M = Vfa(l + k 2 Jk 2 ) 1/2 , 

DiiEMAM = V3tt, DiiM,k z iM = V3Tck z /k. (21) 



Inserting Eq. |21J) into Eq. ( 20 ) and using the leading 
term in the T E lm expansion, we find the asymptotic en- 
ergy between a perfect metal, plasma and Drude cylinder 
and a spherical particle or an atom. Using Eq. ( 171 and 
the T-matrix expansions up to re 3 , see Eqs. ( 10 )- ( 14 ) , we 



obtain the general expression for the asymptotic energy 
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he 
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dn / dk,T EE 



x X (p 2 K 2 {pd) + (2A re 2 + k 2 z )K 2 (pd)) , (22) 

with x — R-s f° r the spherical particle and x = c( E (k) for 
the isotropic atom. Moreover, Aq is proportional to T S A {, 



see Eq. (20), with Aq = 1/4 for the perfect metal parti- 
cle, A = 2| (27ri? s /A^)/(4ii(27riVA^)) for the plasma 
particle and A = for the atom and the Drude particle. 
(17) The latter is due to the fact that T S A [ is set to zero in 



Eq. (|22|. T S A { for the atom is indeed zero and for the 



Drude particle scales with re , see Eq. (12). Therefore, 



for the Drude particle the asymptotic energy given by 
Eq. (22) needs a correction because of the re 4 terms in 
Eqs. (1121) and (fill), which is 
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Inserting Eq. ([9| into Eq. ( 22 ) and using the polar coor- 
dinates, re — > pcos(9)/d and k z — ¥ psm{6) / d, we find 
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ir/2 
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i + G{ P ,e) 
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where for the perfect metal cylinder C(p, 9) = 0, 
for the plasma cylinder C(p,9) « £cos 2 (#) with 
£ = Xl/(2n 2 R 2 \n(2d/R c )), and for the Drude cylinder 
C(p,8) = £'cos(0)/p with £' = A CT d/(4vr 2 i?2 ln(2d/i? c )). 

The correction S£ds for the Drude particle in polar 
coordinates (p, 9) reads 
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Now we use Eqs. (24) and (25) for different material 



properties and calculate the Casimir interaction for var- 
ious limiting cases. 



C. The Casimir interaction between a wire and a 
spherical particle 

Below we present the large separation asymptotic en- 
ergies between a metallic spherical particle and a metallic 
wire. 



1. The perfect metal wire 

For a perfect metal wire and a perfect metal, plasma 
or Drude particle, the energy integral in Eq. (24) results 
into 
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It is important to note that Eq. ( 26 1 depends on the 



material properties of the spherical particle through the 
quantity A . For a perfect metal wire and a plasma par- 
ticle, in the limiting case of small plasma wavelengths, 
A' p -C R s , we reproduce the perfect conductivity form 
with Ao « 1/4. In the limit of large plasma wavelengths, 
A'p > i? s , we obtain A « 7r 2 i? 2 /(15A'p). Since A < 1, 
the plasma wavelength of the spherical particle does not 
have a significant contribution to the asymptotic energy. 
For the Drude particle, using Eq. (25), the correction 



to the asymptotic energy given by Eq. ( 26 1 reads 



S£ps 
lie 



32 d 5 ln(2d/i? c ) V 2A^ 64 ir 



n 2 R 2 



63 A' 



(27) 



and scales with d~ 5 . Depending on the terms in the 
brackets, this correction can have a significant contri- 
bution to the asymptotic energy. In the limit of high 
conductivity, A^ <C R s , the first term in the brackets 
dominates over the second one. For this specific case, if 
X' a d <C R 2 the correction becomes even larger than the 
asymptotic energy itself, i.e. <5£ds 3> £■ 



Note that in Eq. (27 1, the second term in the brackets 
dominates only at low conductivity limit 3> R s in 
which the spherical particle is considered to be a very 
poor conductor. In general, the second term does not 
have a noticeable contribution to the asymptotic energy 
for good conductors such as copper and gold. 



2. The plasma wire 

For a plasma wire and a plasma or a perfect metal 
particle, the polar and radial integrals in Eq. (24) can 
easily be performed, 



Rt 
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he ird 4 \n(2d/R c ) 



f(0 



(28) 



with 

/(0 = ~(2A -1) i 1 (i + O-^ + ^l+O"' > (29) 

where £ is given 
A 2 /(27r 2 i? 2 ln(2d/i? c )). 

length limit, X p /R c -C y/]n(2d/R c ), we reproduce the 
perfect metal wire results given by Eqs. (26) and (27). 

In the opposite limit X p /R c 3> y // \n(2d/R c ), the 
asymptotic energy reads 



below Eq. (|24j), f = 
In the small plasma wave- 



in the case of a perfect particle and a plasma wire. In the 
opposite limit A p > R s , we have A « 7r 2 i? 2 /(15A'p) <C 1 
and the plasma wavelength of the particle does not have 
a significant effect on the asymptotic energy, and the 
asymptotic energy is mainly dominated by the material 
properties of the plasma wire. 

For a Drude particle, the asymptotic energy is given by 
Eq. (30) together with a correction obtained by Eq. (25). 



The correction reads 
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(31) 

Since for good conductors A^ -C R s , similar to the case 
of a perfect metal wire, the second term dominates over 
the first one. 



3. The Drude wire 

For a Drude wire and a perfect metal or plasma parti- 
cle, in the limit f < 1 or d 2 /R 2 < d/X a , using Eqs. d24h 



and (25 1, we reproduce the perfect metal wire results, see 
Eqs. ([26]) and 



In the opposite limit, £' ^> 1 or equivalently d / R 2 



d/X a and d 3> X^/Xa the integrations in Eq. (24) result 
into 
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(32) 



For a perfect metal particle, Ao = 1/4, and the energy 
given by Eq. (32) is always attractive since X a d » R 2 . 



For the plasma particle, in the limit of small plasma 
wavelengtha X' p -C R s , A ~ 1 and the particle behaves 
like a perfect metal. In the opposite limit, X' p >• R s , 

as previously seen, we find Ar ~ ,>J 



case the second term in Eq. (32) dominates, which means 



o w 7r 2 i?f/(15A'p). In this 



that the material properties of the plamsa particle does 
not have a significant effect on the asymptotic Casimir 
energy. 



Under the same condition £' 3> 1, using Eq. (25), the 
correction to Eq. ( 32 ) reads 
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As discussed above, for metallic particles, the second 



term in Eq. ( 33 ) is much smaller than the first one. 
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For a perfect metal spherical particle, Ao = 1/4, and 
the second term can be neglected as we are at the large 
plasma wavelength regime. 

For a plasma spherical particle with plasma wavelength 
A p , if A p <C -R s , A ~ 1/4 and the energy is the same as 



D. The Casimir energy between a wire and an 
Atom 

We calculate the Casimir energy between a wire and 
an atom in both retarded and non-retarded limits. To 
find the asymptotic energies at large separations, we use 
Eq. |24| with A = 0: 



G 



1. The retarded limit 

In the retarded limit, d 3> dio, we find the Casimir 
energy between a perfect metal wire and an atom as 
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37rd 4 ln(2rf/i? c ) 



(34) 



Equation ( 34 ) is in complete agreement with the results 
in Refs. [20- 23 . It is important to note that even though 
the asymptotic energies for an atoms given in Eq. (34) 



and for a perfect metal particle given in Eq. ( 26 ) have the 



same scaling behavior, the numerical coefficients do not 
match; l/(37r) for the atom and 5/ (127r) for the spherical 
particle. This discrepancy is due to the lack of magnetic 
polarizability in the isotropic atoms. 

For the plasma wire and an atom in the limit £ 3> 1 
or equivalently X p /R c 3> ^/ln(2rf/i? c ) and d >• dw, the 
asymptotic energy reads 
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(35) 



which is in agreement with the atom-plasma wire result 
in Ref. [20] ■ In the limit X p /R c < y/ln(2d/R c ) and d > 
dio, we reproduce the perfect metal wire-atom interactin 
energy given in Eq. ( 34 ) . 



For the Drude wire and an atom, in the region of inter- 
mediate distances £' 3> 1 or equivalently d/R c S> R c /X a , 
in the retarded limit d 3> d\o, the asymptotic energy 
reads 
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Equation ( 36 1 is in agreement with Ref. [2UJ • I n the op- 
posite limit, £' <C 1 or d/R c <C R c /X a and d dio, once 
again we reproduce the perfect metal wire-atom asymp- 
totic interaction energy, see Eq. (34). 



2. The non-retarded limit 



In the non-retarded limit, d <C g?ioj using Eq. (24) with 



A = 0, we perform the angular integral for the perfect 
metal, plasma and Drude wires. 

For the perfect metal wire, the radial integral can easily 
be obtained. Expanding the result of the integral for 
d <C dio yields 
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16 d 3 ln(2rf/i? c ) 



(37) 



For the plasma wire, the radial integral over p in 
Eq. ( 24 ) cannot readily be performed in the non- retarded 
limit. Therefore, we expand the integrand for d/d±o <C 1. 
In the limit d <C R c dio/X p the expansion of the inte- 
grand does not depend on the material properties up to 
the leading order. Therefore, performing the integration 



over p results into the perfect metal wire-atom interac- 
tion energy given in Eq. ( 37 ) . 



Similar to the plasma wire, the radial integration in 



Eq. (24) is not easily calculable for the Drude wire. Anal- 



ogously, we expand the integrand for d/d\o <C 1 and then 
perform the integral over p. In the limit d -C R c U d\ojX ai 
we find Eq. ( 37 1 for a perfect metal wire and an atom. 



This is due to the fact that the material properties of 
the Drude wires does not play any role ate intermediate 
separations, see Refs. [TTl 112) . 



E. Universality 

In previous sections, we have derived the asymptotic 
energies between a metallic wire and a metallic spherical 
particle for different dielectric properties, described by 
the Drude, plasma or perfect metal models. We have 
found that in all cases, the Casimir energy depends on 
the material properties of the spherical particle. This is 
due to the fact that the material property of a sphere has 
a significant contribution to its T- matrix |15) . 

In contrast, for parallel metallic wires and a wire-plate 
geometry, at intermediate distances, the asymptotic en- 
ergy does not depend on the material properties of the 
objects and is universal [TTl 112] , 

Although in a wire-sphere system the Casimir interac- 
tion depends on the material properties of the particle at 
all separations, the signatures of the universal behavior 
of the metallic wire is still traceable at asymptotic sep- 
arations. For a plasma wire at intermediate distances, 
d/R c -C exp(Ap/i? 2 ), the Casimir interaction depends 
both on the material properties of the wire and the par- 
ticle. For larger separations, d/R c >• exp(X 2 /R' 2 ), the 
interaction is independent of the material properties of 
the plasma wire, while it still depends on the material 
properties of the particle, see Fig[T^,. 

For a Drude wire and a metallic particle, at larger sep- 
arations, d 2 /R 2 >• d/X a , the interaction depends both 
on the material properties of the Drude wire and that of 
the spherical particle, see Fig[T]3. At intermediate dis- 
tances d 2 /R 2 -C d/X a , while material properties of the 
sphere has a significant role in the Casimir energy, it does 
not depend on the material properties of the wire. Note 
that the casimir interaction is also independent of the 
material properties for two parallel Drude wires [111 112) . 



IV. INTERMEDIATE-SEPARATION REGIME: 
NUMERICAL CALCULATIONS 

We use Eq. pi), to numerically calculate the Casimir 
energy. The numerical algorithm consists of three major 
parts: (i) constructing the matrix N from Eq. (ii) 
computing the determinant of 1 — N for specific imagi- 
nary frequencies k and (iii) integrating over k. The ma- 
trix N consists of blocks which are associated with the 
quantum numbers t and £'. We truncate £ and I' at a 
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51 

T3 




Ap/i? c 



ft? 

t3 



-S?!r§(5Ao + 41n(A^/^)) 




3tt d 4 ln(2d/H c ) 



FIG. 1: Interaction between a wire and a particle and be- 
tween a wire and an atom in the retarded limit. The for- 
mulas describe the rescaled interaction energies, £/(hc). The 
parameter xo = a o for the atom and Xo = Rs f° r the par- 
ticle, (a) Interaction between a plasma wire and a metallic 
spherical particle. The asymptotic results corresponds to the 
regimes sufficiently far from the curve ln(d/i? c ) ~ (X p /R c ) 2 
and for d/R c , X p /R c 3> 1. (b) Interaction between a Drude 
wire and a metallic particle. The separating curve is given, 
up to logarithmic corrections, by d/Rc ~ *J d/X a . The en- 
ergy expressions hold for d/Rc, d/X a 3> 1 and d 3> \\j\ a - 
Different regimes can be reached depending on the relative 
size of length scales: Arrow (1) corresponds to an increasing 
distance d which ultimately leads to a strictly non-universal 
interaction. Arrow (2) indicates an overall increase of the ge- 
ometry (i.e., d/R c fixed) with constant conductivity leading 
to a " universal " interaction. 



finite partial wave number £ majc such that the result for 
the energy changes by less than a factor of 1.0001 upon 
increasing £ max by 1. Since t, £' < £ max , the N-matrix has 
^max blocks N u' ■ The block N^' consists of the elements 
N/ m f' m i which form 2x2 blocks, with m — —I, ...,(. 



and m! 



N 



£m,£'m' 



mEAl mEE 



Consequently the size of the block N«< is {At + 2) x 
{AC + 2), implying that the off-diagonal blocks {£ ^ £') 
are not square matrices. Furthermore, N«/ blocks are 
not diagonal since symmetry along the axis parallel to 
the wire's axis is broken by the particle. 

To construct the matrix N, 4£ 2 nax (i' max + 2) 2 integrals 
over k z have to be evaluated for each n. This makes 
the numerical computations for closer separations quite 
expensive. For example at separation d/R — 2.6, the en- 
ergy converges with £ ma , x = 12, corresponding to a matrix 
of size 336 with 112896 fc z -integrals for just a single value 
of K. 

Figure [2] illustrates our numerical results for a metallic 
cylinder and a metallic spherical particle both with the 
radius R. The plots show the Casimir energy normal- 
ized to the energies for the perfect metal wire-particle 
configuration as a function of the surface-to-surface dis- 
tance h = d — 2R. For the numerical calculations, wc 
used Xp/R = 0.05 and 0.5 with X p /\ a = 27.4, corre- 
sponding to the parameters for gold with A p = 137 nm 
and A CT « 5 nm [53]. Figure [2] shows the dependence of 
the Casimir energy on the material properties of the wire 
and particle. Figure [2]ja) depicts the interaction energy 
between the perfect metal wire with the metallic parti- 
cle for plasma (opensymbols) and Drude (filled symbols) 
models with X p /R = 0.05 (squares) and 0.5 (circles) and 
A p /Act = 27.4. The dashed lines show the asymptotic en- 
ergies given by Eqs. ( 26 ) and ( 27 1 . As shown in the figure, 



there is a very good agreement between the asymptotics 
and the numerical results. The numerical results confirm 
that at very large separations, only the material proper- 
ties of the particle contribute to the Casimir energy. We 
note that in the range of distances considered here and 
for the parameters of gold, one has to consider the en- 
ergy given in Eq. ( 26 ) in addition to the correction term 
presented in Eq. (27 1. 



Figure |2jb) shows the interaction between a plasma 
wire and a perfect metal (open triangles), plasma (open 
circles) and Drude particles (filled circles) with the 
plasma wavelength X p /R = 0.5 and A p /A CT = 27.4. 
The dashed lines show the asymptotic energies given by 
Eqs. (|30|) and (pTl). This figure also shows the good 



agreement between the asymptotics derived in the pre- 
vious section and our numerical results. Again, here it 
is important to include the correction given by Eq. (31 1. 



Note that for the perfect metal particle and the plasma 
cylinder (open triangles) the energy ratio approaches 1 
at large separations. This is due to the fact that in this 
regime the material properties of the plasma wire do not 
contribute to the Casimir interaction, see Fig. [lja) . 

Figure^c) shows the interaction between a Drude wire 
and a perfect metal (open triangles) , a plasma (open cir- 
cles) and a Drude spherical particles (filled circles) with 
Xp/R = 0.5 and Ap/A CT = 27.4. The dashed lines are ob- 



tained by computing the integrals in Eqs. (|24|) and (251 



Since for 10 2 < d < 10 4 we have 10~ 2 < f < lO" 1 , corre- 
sponding to the crossover regime, the asymptotic energies 
of Eqs. (32 1 and (33) are not applicable, see Fig. [If b). 
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FIG. 2: Ratio of the numerically computed energy E for re- 
alistic metals to the perfect metal energy E pal { versus h/R 
for a wire and a particle with R c = R E — R. Each fig- 
ure shows the numerics related to a wire with specific ma- 
terial properties: (a) for a perfect metal wire; (b) a plasma 
wire with X p /R = 0.5 and (c) a Drude wire with the same 
value for \ P /R and X a = A p /27.4. The ratio of energies for 
particles with different material properties are distinguished 
by different symbols: Open triangles for perfect metal par- 
ticle and open circles and open squares for plasma particles 
with Xp/R = 0.5, 0.05, respectively. Filled circles and filled 
squares denote Drude particles with X p /R — 0.5, 0.05, re- 
spectively, and X a = Ap/27.4. 
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FIG. 3: The geometry of a wire of radius R c and a sphere 
of radius R s at a distance ho (center-to-center distance d = 
h + R C + R s ). 



V. SHORT-SEPARATION REGIME: 
PROXIMITY FORCE APPROXIMATION 

In this section, using the Proximity Force Approxima- 
tion (PFA) [25], we calculate the Casimir interaction at 
short separations h <C R c , Rs- This method gives the in- 
teraction as an integral of the energies between parallel 
surface segments, 



£pfa = / dAE p i ate (h) . 



(38) 



where -Spiatc is the Casimir energy per unit area between 
two parallel plates and h is the surface-to-surface dis- 
tance. Figure [3] illustrates the distance between two sur- 
face elements. According to Fig. [3] the distance h is 

h = h + R S {1 - cos(0 s )) + R c (l - cos(0 c )) , (39) 

with ho = d—R c — R s the distance of the closest approach 
between the cylinder and the sphere. One can write <f> c 
in terms of 9 S and 4> s , 

sin(0 c ) = ^ sin(0 s ) sin(0 s ) . (40) 

At short separations, the surface elements of the sphere 
and cylinder in which 9 S <C 1 and <f) c <C 1 respectively, 
contribute most to the interaction. Therefore, the dis- 
tance h can be approximated by 



^o + ^(l + f sin 2 (0 s )). (41) 



Inserting Eq. (41) into Eq. (38) and performing a simple 

e PFA energy 

dHE platc (H). (42) 



change of variable, we obtain the PFA energy 

OPFA — 



For the case of perfect metal surfaces with E p i at0 (H)/A = 
~Tr 2 hc/(720H 3 ), we find 

£ppa = . (43) 
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FIG. 4: Ratio of the numerical results for the Casimir en- 
ergy shown in Fig. [2] and the PFA energy for perfect metals 
(squares) and for plasma model with X p /R = 0.5 (circles). 
The ratio is shown as a function of surface-to-surface distance 
h. Similar to Fig. [2j the wire and particle have equal radiuii 
Rs — Rc — R- 



For Drude and plasma objects, we use the Lifshitz for- 
mula US, 



E p \ ane (h) _ he 



pdp 



A (2tt) 2 
x In [(1 - rf r^e- 2Kph ){l - rf r^e~ 2Kph )] , (44) 



with p a dimensionless variable, and r 1 " ylJ> the 

Fresnel coefficients of the surface elements for the cylin- 
der and sphere, respectively. The Fresnel coefficients of 
an object a are given by 



M (E) 



p. a {icn) - V 1 " 


f (n 2 a (ic K ) 


- l)/p 2 


Mo(*CK) + y/l - 


\- (n 2 a (icn) 


- l)/p 2 


e a (icn) 


- (n 2 a (icK) 


- w 


e a {icn) + V 1 "t 


- (n 2 a (icK) 


- l)/p 2 



the 



refraction index, n a (icn) 



(45) 



r^(icK,p) -- 

r a(iCK,p) -- 

where n a is 
\J e(icn) /i(icn) . 

The PFA energy is obtained using Eqs. (|42|) and ( |44[ ) 
together with the dielectric function of Eq. Q. Figure 
[4] shows the Casimir energy for a perfect metal (open 
squares) and plasma model X p /R = 0.5 (open circles) 
normalized to the PFA energy. The energies associated 
with the Drude model are not shown since they collapse 
on the data for the plasma model at short separations. 
Our data show that as the distance between the sphere 
and the cylinder decreases, the values of the energies be- 
come closer to the PFA ones. Note that at short sepa- 
rations the energy converges with larger values of £ max - 
Since the size of the matrix N increases quadratically 
with £ max , the numerical calculation of the Casimir en- 
ergy becomes extremely costly at short separations. In 
this work the interaction is calculated up to h/R — 0.4. 



VI. SUMMARY AND CONCLUSIONS 

In summary, we have studied the Casimir energy be- 
tween a cylindrical wire and a spherical particle. For 
large separations, we have derived the asymptotic ener- 
gies for the Drude, plasma and perfect metal models. In 
addition, we have calculated the Casimir interaction be- 
tween a metallic wire and an isotropic atom. Our results 
for the wire-atom system is in complete agreement with 
previous results obtained through a different method [2DT 

• 

Furthermore, we have computed the Casimir interac- 
tion between a spherical particle and a wire. Such com- 
putations are quite demanding due to lack of spherical 
symmetry. Our numerical results perfectly match the 
asymptotic energies. 

For short separations, we obtained the energy using the 
Proximity Force Approximation (PFA) and compared it 
with our numerical data. This comparison indicates that 
as the distance between the wire and particle decreases, 
the numerical results for the Casimir energy becomes 
closer to the PFA one. It is noteworthy that depending 
on the separation, the material properties of the metallic 
wire may not play a role in the interaction energy, similar 
to the parallel wires and wire-plate systems [TTJ [T2] • 

In a cylinder-sphere system, we do not observe a uni- 
versal behavior as we have previously obtained for par- 
allel wires and a wire-plate geometries because of the 
physical properties of the spherical particle. However, 
one can still have "universal" regimes in which the inter- 
action does not depend on the material properties of the 
metallic wire. 

In case of the plasma wire with the plasma wave- 
length A p and radius i? c , at sufficiently large separations, 
d/R c ^> exp(A 2 /i? 2 ), the material properties of the wire 
does not play any role in the asymptotic interactions be- 
tween the wire and a particle or an atom. 
In contrast for the Drude wire with conductivity a and 
the characteristic length A CT = 27rc/er, at large separations 
d 2 /R 2 ^> d/X a , the asymptotic energy depends on the 
the material properties of the wire. Quite interestingly, 
in the opposite limit, d 2 /R 2 <C d/\ a and d ^> X 2 /X a , the 
asymptotic interaction becomes independent of the ma- 
terial properties of the Drude wire. The specific behavior 
of the Drude wires has been explained in Refs. [HJ[T2] in 
terms of large scale charge fluctuations. 

At the end we emphasize that since simple generic ge- 
ometries appear in many nano- and micrometer-sized 
systems, the knowledge of the interaction between a 
metallic sphere and cylinder could be important for an 
efficient design of low-dimensional structures. 
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No. S-000354 (RZ and TE). The T-matrix of a sphere in spherical vector wave ba- 

sis is diagonal in the quantum numbers I, m and the 
electromagnetic polarizations E and M 

Appendix A: T-matrices 

In this appendix, for completeness, we present the T- 
matrices of a cylinder and a sphere, from Refs. |1U 115) . 



n t)I 1+ x(kR) 


h+± { n K>R) + 2nnRI' l+ 1 (nnR) 


— nl l+ i (nnR) 




-2kRI'AkR) 

'2 




2 r,K l+i ( K R) 


I l+ i (nnR) + 2nnRI' l+ s (nnR) 


- nI l+ i(nnR) 




f 2kRK' i (kR) 

t+ 2 



(Al) 



with n = yf e(in)fj,(iK) and n = J e(in) / JjJJk). The T- 
matrix elements for E-multipoles, T E 



Imi 



are obtained 
from Eq. ( Al I by interchanging e and p. 

The T-matrix elements of a cylinder with dielectric 
response t(in) and magnetic permeability p(in) are given 
by 



where Ce = —Cm = 1/(2ttL) is the normalization coef- 
ficient with L the overall length of the cylinder. 
Inserting Eq. (|B2|) into Eq. (|Blh yields 



\E) 



\ rj k'm'P' 



E 13 



out \ thPP' 
k z mPI s.k z m,k'm' i 



(B3) 



rpEE __ 



InipR) A 2 A 3 + if 2 



T, 



EM 

k z n 



K n (pR) AxA 2 
K 



K 2 ' 



with X = {nk z /(^/IpR 2 n)) (l/p 12 - 1/p 2 ) and 



4(p'i?) 1 ^(pi?) 



p'RI n (p'R) e pRK n (pR) 



(A2) 
(A3) 

(A4) 



with 



T a P £ m , Km , = (-l)C P (E^ mP \T s \E^ m , pi ) , (B4) 

the T-matrix of a sphere in cylindrical basis. Now we 
expand the cylindrical basis wave functions in terms of 
the spherical basis waves, 



\El 



2_j D traQ,k,m,p\EY^Q 



(B5) 



Note that A 2 can be obtained from Eq. (A4) by inter- 
changing e with (i, and A3 can be fou nd by replacing K' n 
with I' n and K n with 7„ in Eq. ( |A4| ). Moreover, T^f 
can be obtained by replacing e with and considering 



rpME 
± k z n 



_ _rpEM 
~ ± k z n ■ 



Appendix B: T-matrix of a sphere in cylindrical 
basis 



The electromagnetic field far enough outside a sphere 
can be written in terms of the regular wave function 
with the electromagnetic polarization P, \E/ m p(icK)) j 
the free electromagnetic Green's function, G (zck), and 
the scattering operator of the sphere, T s (icn) [7], 



\E) 



I B reg \ 

\ rj k z mP> 



TJE 



k z mPl 



(Bl) 



where icn arguments are dropped for brevity. 

The expansion of the free Green's function in terms of 

the regular and outgoing wave functions is given by 



E 

k z mP 



C P \EZ mP )(ElJ 



mP \ 1 



(B2) 



where Q is the electromagnetic polarization, I is the 
quantum number related to the spherical wave functions 
and the coefficients Di m Q^ z mP are the elements of the 
conversion matrix from the cylindrical to spherical basis, 
see Appendix [C] for the detailed description. Since the 
azimuthal dependence of the wave functions in both sides 
of Eq. (B5) are the same, the sum runs on the quantum 



number I and polarization Q. 



Inserting Eq. (B5) into Eq. (B4), we obtain 



T. 



with 



pp' 



z m,k'm 



/ 1 C" ( v\ k z mP,tmC 
>Q> 



£Q,£'Q 



xt; 



QQ' 



Iml'm 1 Dt'm'Qi \k'm'P' i 



lml = {-l)C' Q {.)(E\Z Q \^\E\Z' Q , 



(B6) 



(B7) 



the T-matrix of the sphere in the spherical basis, see Ap- 
pendix |Aj and C' m (k) = —C' e (k) = K the normalization 
coefficients of the Green's function expansion in spheri- 
cal basis. The ratio of the normalization coefficients in 



11 



Eq. p6]) is (1 - 28 p ,q)/(2ttkL). Since the T-matrix of 
the sphere in spherical basis is diagonal in /, m and po- 



larization, Eq. ( B6 ) is simplified to 
1 



T 



pp i 



s,k z m,k> z m 2nKL 
t 



J2 E(!- 2 ^q) 

max(l,|m|) Q 



xD LrnP,£ m Q T 2lrn D imQ,k' z mP' ■ (B8) 



Appendix C: Conversion Matrix T)tm,k z m 

The coefficients of the expansion of the cylindrical vec- 
tor waves in terms of the spherical ones determine the 
elements of the conversion matrix Di m .k z m- These coef- 
ficients are known and have already been calculated [271 - 
129] . In this appendix we make the previously derived 
coefficients consistent with the Wick-rotated vector wave 
bases introduced in Ref. [7 . The expansion of cylindri- 
cal vector waves (m m ^, n m \) in terms of spherical vector 
waves (m£ m ,ng m ) is given by [2"5] 

oo 

l—m 

oo 

n m A = 2^ Aem\nim + Be m \me m , (CI) 



where A 2 = k 2 — h 2 and 



.4 



2^+1 



m)\ t _ 



e(i+i) (e + m)\ 



m+1 



fcsin(a) — P» m (cos(a)) , 
da 



with cos(a) = h/k. Note that m and n are the vector 
wave functions in Euclidean space. 

Taking into account h = k Zl k = in and A = 
«a/k 2 + k 2 , the vector wave bases (M rcg , N rog ) defined in 
Ref. [7] are related to the bases (m, n) defined in Ref. [28] 
by the relations 



n mA = * m - 1 v / ^Tfc 2 N^ s m ( K ,x) 



(C3) 



and 



1 V {2l + l)(£-m)\ lm 



mL e (m), 



Plugging Eq. \G2\ (after a Wick rotation [k — > in)), 
Eq. (C3) and Eq. (j^J) into Eq. flC!] ), we obtain 



l—m 



+^mE,fc z mMN2(^,x) , 

oo 



+AmE,fc 8 m£N^(K, x) 



(C5) 



where the conversion matrix elements read 



D 



imM,k z mM 



= (-!)< 



1.2 \ 2 



tmE,k z mM 



Utt(2£ 


■f l)(£-m)\ 


/ + 


l){£ + m)\ 


2 p'T(- 


-ik z /n) , 


Utt(2£ 


■f l)(£-m)\ 


/ + 


l)(£ + m)\ 



hn f 1 + 5) V(-'*«/k). 



DimM,k z mE — — DtmE,k z mM , 
Di m E,k-mE = DtmM,k z mM ■ 



Since it is difficult to deal with the Legendre functions 
with complex arguments, we use the Rodrigues represen- 
tation of Legendre polynomials and find 



pp(-ix) = (~iy +m fr(x) 

where the real function /, m is given by 



(C7) 



(C8) 



Using Eq. (C7), we can write the conversion matrix 
k m m terms of a modified matrix D^ m /- m) 



D 



£m,k z m 



(-ir m H) 



m / -\^+m— 1 



D 



£m,k z m •> 



(C9) 



with 



lmM,k z mM 



D 



ImE.k-mM 



Uix(2£ 




— m)\ 


' £(£ + 




m)\ 


1 ^ z 

K 2 t 


1 

\ 2 

\ -fral 

) h 


(k z /n) 


Un{2£ 




— m)\ 



x m I 1 + 



' + l)(£ + m)! 

1 

A, 



(CIO) 



fT{k z /n). 



DtmM,k z mE — —DlmE,k z mM , 
DlmE.k z mE = DtmM,k,mM ■ 
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